25 research outputs found
Progressive construction of a parametric reduced-order model for PDE-constrained optimization
An adaptive approach to using reduced-order models as surrogates in
PDE-constrained optimization is introduced that breaks the traditional
offline-online framework of model order reduction. A sequence of optimization
problems constrained by a given Reduced-Order Model (ROM) is defined with the
goal of converging to the solution of a given PDE-constrained optimization
problem. For each reduced optimization problem, the constraining ROM is trained
from sampling the High-Dimensional Model (HDM) at the solution of some of the
previous problems in the sequence. The reduced optimization problems are
equipped with a nonlinear trust-region based on a residual error indicator to
keep the optimization trajectory in a region of the parameter space where the
ROM is accurate. A technique for incorporating sensitivities into a
Reduced-Order Basis (ROB) is also presented, along with a methodology for
computing sensitivities of the reduced-order model that minimizes the distance
to the corresponding HDM sensitivity, in a suitable norm. The proposed reduced
optimization framework is applied to subsonic aerodynamic shape optimization
and shown to reduce the number of queries to the HDM by a factor of 4-5,
compared to the optimization problem solved using only the HDM, with errors in
the optimal solution far less than 0.1%
An adaptive, training-free reduced-order model for convection-dominated problems based on hybrid snapshots
The vast majority of reduced-order models (ROMs) first obtain a low
dimensional representation of the problem from high-dimensional model (HDM)
training data which is afterwards used to obtain a system of reduced
complexity. Unfortunately, convection-dominated problems generally have a
slowly decaying Kolmogorov n-width, which makes obtaining an accurate ROM built
solely from training data very challenging. The accuracy of a ROM can be
improved through enrichment with HDM solutions; however, due to the large
computational expense of HDM evaluations for complex problems, they can only be
used parsimoniously to obtain relevant computational savings. In this work, we
exploit the local spatial and temporal coherence often exhibited by these
problems to derive an accurate, cost-efficient approach that repeatedly
combines HDM and ROM evaluations without a separate training phase. Our
approach obtains solutions at a given time step by either fully solving the HDM
or by combining partial HDM and ROM solves. A dynamic sampling procedure
identifies regions that require the HDM solution for global accuracy and the
reminder of the flow is reconstructed using the ROM. Moreover, solutions
combining both HDM and ROM solves use spatial filtering to eliminate potential
spurious oscillations that may develop. We test the proposed method on inviscid
compressible flow problems and demonstrate speedups up to an order of
magnitude.Comment: 29 pages, 13 figure
A space-time high-order implicit shock tracking method for shock-dominated unsteady flows
High-order implicit shock tracking (fitting) is a class of high-order,
optimization-based numerical methods to approximate solutions of conservation
laws with non-smooth features by aligning elements of the computational mesh
with non-smooth features. This ensures the non-smooth features are perfectly
represented by inter-element jumps and high-order basis functions approximate
smooth regions of the solution without nonlinear stabilization, which leads to
accurate approximations on traditionally coarse meshes. In this work, we extend
implicit shock tracking to time-dependent problems using a slab-based
space-time approach. This is achieved by reformulating a time-dependent
conservation law as a steady conservation law in one higher dimension and
applying existing implicit shock tracking techniques. To avoid computations
over the entire time domain and unstructured mesh generation in higher
dimensions, we introduce a general procedure to generate conforming,
simplex-only meshes of space-time slabs in such a way that preserves features
(e.g., curved elements, refinement regions) from previous time slabs. The use
of space-time slabs also simplifies the shock tracking problem by reducing
temporal complexity. Several practical adaptations of the implicit shock
tracking solvers are developed for the space-time setting including 1) a
self-adjusting temporal boundary, 2) nondimensionalization of a space-time
slab, 3) adaptive mesh refinement, and 4) shock boundary conditions, which lead
to accurate solutions on coarse space-time grids, even for problem with complex
flow features such as curved shocks, shock formation, shock-shock and
shock-boundary interaction, and triple points.Comment: 35 pages, 20 figure
Accelerated solutions of convection-dominated partial differential equations using implicit feature tracking and empirical quadrature
This work introduces an empirical quadrature-based hyperreduction procedure
and greedy training algorithm to effectively reduce the computational cost of
solving convection-dominated problems with limited training. The proposed
approach circumvents the slowly decaying -width limitation of linear model
reduction techniques applied to convection-dominated problems by using a
nonlinear approximation manifold systematically defined by composing a
low-dimensional affine space with bijections of the underlying domain. The
reduced-order model is defined as the solution of a residual minimization
problem over the nonlinear manifold. An online-efficient method is obtained by
using empirical quadrature to approximate the optimality system such that it
can be solved with mesh-independent operations. The proposed reduced-order
model is trained using a greedy procedure to systematically sample the
parameter domain. The effectiveness of the proposed approach is demonstrated on
two shock-dominated computational fluid dynamics benchmarks.Comment: 24 pages, 8 figures, 2 table